A new variable in flow analysis 
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We have used a simple spectrum distribution which was derived from a hydrodynamical equation 
[12] to fit the data of the STAR group. It is found that it can fit the V2 of STAR group very well. 
We have found that V2 is sensitive to both the effective temperature of particles and the expanding 
velocity. We have suggested a new variable z to be used in the flow analysis. This new variable will 
measure the correlation of particles momentum components. We have also shown that one of the x 
or y direction in the reaction plane is the direction which has the largest variance. 



Recently, the study of elliptic flow has attracted atten- 
tion of both theoreticians and experimentalists. [1-8]. It 
was argued that an anisotropic distribution of final state 
particles with respect to the reaction plane can be used to 
reflect the strong re-interaction among quarks and gluons 
in the initial state. 

The calculation of t>2, the anisotropic distribution of 
final state particles is always done in the following two 
ways:(l): we determine the reaction plane first, then we 
calculate the average of (cos(2(</> — Here <j)R is 

the azimuthal angle of the reaction plane and 4> is the 
azimuthal angle of particles in the Lab frame. (2): us- 
ing the pairwise azimuthal correlations [11] to calculate 
the i>2- The second method has an advantage over the 
first method that no reaction plane is needed to be de- 
termined. But a prior knowledge on the distribution of 
P(4>) is needed to be known for the second method. The 
biggest problem for the first method is to determine the 
reaction plane. There are several methods are used in 
data analysis. In RHIC, the so called second harmonic 
event plane is widely used in the data analyzes. 

However, it seems to the authors that we still not clear 
what we really have measured in the experiment. To dis- 
cuss this question, we will use the spectrum distribution 
in Ref. [12]. This spectrum distribution is a solution for 
a hydrodynamical equation and it can be expressed as 
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Here v = (v x ,v y ) is the expanding velocity of particles. 
T x and T y are effective temperatures of x component and 
y component of particles momentum [12]. The above pic- 
ture can be understood in the following way: A A collision 
will form a QGP or a dense hadron phase at the initial 
state. Due to the initial asymmetry collisions, it is ex- 
pected that the effective temperature will be different in 
x and y components. After some time, the hadron gas 
are formed. This hadron gas will expand in velocity v. 
This distribution is relative to the reaction plane. 

Using this equation, we will calculate v-i using fol- 



lowing three variables. They are p, p = p p, and 

z = ( ' Wfrj)- Here and are second mo- 
ments in x and y direction respectively which are defined 

as 
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Then we will calculate vi using the following formula 
f P(u)%4du 
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Here u is one of the three variables in the above formula. 
It is easily checked that 



for all cases. 
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when T x — T y ; but V2{p) is not zero when T x ^ T v . i?2(p) 
is not zero and its value is shown in Fig.l. It is interest- 
ing to notice that flow may be caused by the expansion 
velocity v and the " effective temperature" if we use vari- 
ables p and p. However, there is no flow to be observed 
if we use variable z. 
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v=(0.273,0.2),T x =0.12 GeV,T y =0.12 GeV 
v=(0.273,0.2),T>0.13 GeV,T y =0.12 GeV 
v=(0.2,0.5),T x =0.12 GeV,T y =0.12 GeV 
Star Data from Ref.[2] 
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Its eigenvectors determine the x direction and y direc- 
tion of the reaction plane if we take u = p. However 
transverse momentum method is widely used in the data 
analysis to determine the direction of reaction plane. The 
reaction plane is determined by using the following for- 
mula: 
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FIG. 1. vi vs. p t . 

From Fig.l, we notice that:(l) the differences between 
T x and T y can cause strong flow. (2): The difference 
between v x and v v also cause flow. It seems that the data 
does not support the assumption that T x is different from 
T y . But we need to point out that this results is model 
dependent. One interesting thing is that this model can 
fit the data very well especially in the larger px region. 

For elliptic flow analysis, we are more interested in 
the distribution of particles in momentum space or az- 
imuthal angle space. We would like to see the relative 
differences between x-components and y-components of 
particles momentum. The difference between T x and T y 
might be caused by the initial flow in QGP phase. But 
T x and T y are measurements which could tell us the abso- 
lutely value of particles momentum in different direction. 
Besides this effective temperatures difference we also like 
to find a correlation inside the distribution of P(p x , Py)- 
From the above simple model, we find that z variable is 
a proper variable for this purpose. We would like it can 
be used in the data analysis. 

To understand the meaning of the flow, we will start 
from the transverse momentum sphericity tensor defined 
by 
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where Ui(i = 1,2) is the i-th component of the variable 
u. As above, this u can be one of three variables mention 
above. To be more clearly, we have the following matrix. 
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It is clear that the eigenvalue of this matrix are 
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Here = 1 when particles rapidity is big than zero and 
uji = —1 when particles rapidity is less than zero. The 
basic idea behind this method is that the total transverse 
momentum of the system in the y direction of the reaction 
plane is zero and the symmetry of the collision system in 
the direction of to O + n. Here is the polar angle. Can 
this two methods give the same reaction plane? 

We will give a positive answer to this question under 
the condition that the multiplicity of event is infinity in 
this paragraph. When we determine the two eigenvec- 
tors, then S12 = in this new frames. If we use Eq.(8) 
to determine the reaction plane first, then we can also 
find that Si 2 will be zero in the frame due to the sym- 
metry of u 2 «-> — u 2 . We will show that there is only one 
vector in the x — y planes which lead the value of S\ 2 is 
zero. Suppose that a new axis which has angle relative 
to the reaction plane is 9. It is easily checked that the 
new quantity in this new frame is 



t4(i/)=u£(i/)cos(0) + 
u'Ju) = -u*{v) sin(<9) 



<<»sin(0) 
f u£(i/) cos(0) 



(9) 



Here the superscript R refers that those quantities are in 
the reaction plane frame. u'(V) is the variable in the new 
frame. Then 
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It is clear that S' 12 is zero only when 9 = or 7r or 
S\i — S 22 . For the first case, it means that S' 12 = 
when it is in the reaction plane frame. When Sf^ = S 22 , 
the shape of u T Su is circle in momentum space. The 
eigenvector for this matrix is arbitrary. Therefore, we 
have shown that the transverse momentum method and 
sphericity tensor give the same reaction plane when the 
multiplicity is huge. 

Then a which is used to measure the asymmetry dis- 
tribution of particles in the momentum can be calculated 
in the following ways. 
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Using Eq(7), we have 
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Taking u = p t , a can be expressed as 
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Here (pi) = Pt are the transverse momentum 

square average for the event. M is the multiplicity for 
the event. Then we have 
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This average is taken over all events. If the transverse 
momentum is fixed, then we have 



a{p t ) = v 2 (p t )- 



(15) 



Therefore, we can calculate ^(.Pt) using Eq.(12) directly. 
If we choose u = z = (z\, z 2 ) = y ? B - 2 ), then we 

find that the matrix for z becomes to 
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When particles rapidity window is huge and number of 
particle is infinity for each event, we can write 

Sn = J P(p)p x p x dp x dp y 

S 2 2 = j P(p)PyPydp x dp y 

Sl2 = j P(p)p x Pydp x dpy (17) 



Therefore 
a(z) = ■ 
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Thus a(z) actually measures the correlation between the 
two component of particles. If the particles are emitted 
randomly, a(z) — 0. On the other hand, if particles are 
emitted always in the a particular direction, say p x = p y , 
then a(z) = 1. It is easily checked that 
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Here t = It is interesting to notice that when t = 1, 
then a(p) = a(z). It is easily seen that a(p) > a(z). If 
t = 0, then a(pt) = 1 which corresponds to the case par- 
ticles are emitted only in the y direction. The advantage 
of a(z) over a(p t ) is clear: it depends on a dimensionless 
variable; it measures the correlations between the com- 
ponents of particles momentum. 

If we take u = p t and y = ( -^fy , ^fy ) , then we will 

have similar expressions as above. The only difference is 
that we need to use varible p t and y to take the place of 
variables p t and z in above expressions. 
If a(z) is very small, we have 
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Thus a(p) measure the ratio between the difference of 
the variance in x and y directions and the sum of the 
variance in x and y directions. The detail information of 
a(z) will not be observed. 

One of interesting things about the matrix of variable 
z is that its eigenvector is always along the direction 

(^2' \7f) 0r (^2' 7§)' ^ ue to tne f ac t t nat tne nn itc mul- 
tiplicity in the collisions, the "estimated reaction plane" 
will fluctuate randomly. Finally, we need to point out 
that this estimated reaction plane is quite different from 
the reaction plane estimated using variable p since they 
are in different variable spaces. 

In the experiment, experimentalists normally used u = 
= (cos(0),sin((/>)) in the analysis. The matrix for x 

is 
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When the system transform to eigenvectors frame, we 
have 
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This is the way to calculate the reaction plane azimuthal 
angle (<Pr) in RHIC data analyses. 

We can construct other varibles to dertermine the high- 
order harmonics reaction plane. For example, when n = 

4, we can take u = (^F, = (sin(2</>), cos(2</>)). If 

we calculate the matrix in the corresponding eigenvector 
frame, we have 
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This is the way to calculate the fourth harmonic reac- 
tion plane as metioned in the Ref. [9]. In general this 
estimated reaction plane should be different from the es- 
timated second harmonic reaction plane since they are 
estimated for different variables. The corresponding a 
for this matrix is V4. 

We have shown in the above that when the multiplicity 
is huge, transverse momentum and sphcrecity matrix will 
give the same "estimated" reaction plane for variable u 
as long as they have symmetry u y <-> —u y . However, we 
know that this estimated reaction plane is quite different 
for different variable u. What is the physical meaning of 
the direction estimated by sphericity matrix and trans- 
verse momentum? We will show in the next paragraph 
that this eigenvectors directions actually gives the direc- 
tion where the variance is the maximum. 

Suppose that u is a vector (2x1 matrix) in the lab 
frame and we will try to find a new frame such that the 
new variable u' has the largest variance along a axis in 
the frame. Suppose a T is a 2 x 2 matrix, then the new 
variable u'will be 

t4 = a?\i, 2 = 1,2. (23) 

Here &f = (an, 0*2). Then its variance is 

M 

- a?uu T ai = a^Sa;, i = 1,2. 

(24) 

if we take a constraint that a^aj = 1. Then we can 
construct the following Lagrange multipliers 

L(ai) = a^a; - A(a; r ai = 1). (25) 

Taking derivation with ai, we have 

dL 

— = 2Ssli - 2A &i (26) 
o&i 

Setting this value to zero, we have 

Sai = Aai, (27) 

and 

M 

^Mf<M=A i; i = l,2 (28) 

This tell us that when we choose one of the eigenvector 
as a\, then its maximum variate will be the eigenvalue. 
We will choose the largest eigenvalue. It is also easily to 



prove that if we choose 02 as another eigenvalue. Then 
the S12 will be zero for variable u' . Therefore we have 
shown here that the one of the eigenvectors actually gives 
us the direction where the variance is the largest. For an- 
other eigenvector which give us another direction which 
has a property that S12 will be zero (in other words, 
its component will be uncorrelated with the first com- 
ponent). The largest variance for our case is the largest 
eigenvalue. 

Conclusions: It has been shown that the V2 measured 
in the data reflects the information on the expanding ve- 
locity, effective temeparature in x and y components and 
correlation between different components of the particles 
momentum. However due to the fact that the correla- 
tion between different momentum of particles are small. 
Therefore we redefine a new variable z which can be 
used in the data analysis. We belived that V2(z) is small 
but can show us the information on the correlation be- 
tween different components of particles. We have also 
shown that estimated reaction plane will be different if 
we choose different variables in the data analysis. The 
two directions in the reaction plane can be undersood as 
the biggest variance directions of the data. One of the 
interesting thing is that, the model can fit the data quite 
well. This model suggests that the effective tempara- 
ture for different directions of momentums are almost 
the same. 
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